Late-time tails of wave maps coupled to gravity 
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We consider the late-time asymptotic behavior for solutions of Einstein's equations with the wave 
map matter. Solutions starting from small compactly supported ^-equivariant initial data with t > 1 
are shown to decay as t~^ 2e+2 ^ at future timelike infinity and as u~^ l+1 ^ at future null infinity. 



I. INTRODUCTION 

In this paper we continue our investigations, initiated in [ij, of the precise quantitative description of the late- 
time asymptotic behavior of self-gravitating massless fields. In [l[ we considered the simplest case of a spherically 
symmetric massless scalar field. Using nonlinear perturbation analysis we showed that solutions starting from small 
^ initial data decay as t~ 3 at timelike infinity and as u~ 2 at null infinity. We also derived a simple analytic formula for 
the amplitude of the late-time tail in terms of initial data. 

Here we study the analogous problem for wave maps which are a natural geometric generalization of the wave 
equation for the massless scalar field. This generalization seems interesting because in the so called equivariant case 
the homotopy index I of the map plays the role similar to the multipole index for spherical harmonics. However, 
in contrast to the decomposition of a scalar field into spherical harmonics which makes sense only at the linearized 
level, it is consistent to study nonlinear evolution for the wave map within a fixed equivariance class. In this sense 
^-equivariant self-gravitating wave maps can serve as a poor man's toy-model of non-spherical collapse. The £ — case 
reduces to the spherically symmetric massless scalar field analyzed in [l[ so hereafter we assume that I > 1. We note 
aside that the £ — 1 case has been extensively studied in the past focusing on the critical behavior at the threshold of 
black hole formation 0, H, 0) IE @1 j however, to our knowledge, the late-time behavior of wave maps coupled to gravity 
i i' has not been analyzed before. 

Using the same third-order perturbation method as in [l[ we show here that for small compactly supported initial 
data the late-time tail of the self-gravitating ^-equivariant wave map decays as £-( 2£ + 2 ) at future timelike infinity and 
as at future null infinity. We also compute the amplitude of the tail in terms of initial data. These analytic 

I ' results are verified by the numerical integration of the Einstein-wave map equations. 

ON : 

II. SETUP 

\o ■ 

o : 

, Let U : Ai — > N be a map from a spacetime g a b) mto a Riemannian manifold {M,Gab)- A pair (U,g a b) is 
■ said to be a wave map coupled to gravity if it is a critical point of the action functional 
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s = I (tt^ - ^g ab d a U A d b U B G AB ) dv, ( 1 ) 
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where R is the scalar curvature of the metric g a b, G is Newton's constant, A is the wave map coupling constant, and dv is 
the volume element on (M,g a b)- The field equations derived from ([1]) are the Einstein equations R a b~ \gabR = 8^GT a b 
with the stress-energy tensor 

Tab = A (d a U A d b U B - \g ab {g cd d c U A d d U B )\ G AB , (2) 



and the wave map equation 



O g U A + T A c (U)d a U B d b U c g ah = 0, (3) 



where Tq C are the Christoffel symbols of the target metric Gab and D g is the wave operator associated with the 
metric g a b- As a target manifold we take the three-sphere with the round metric in polar coordinates U A = (F, ft) 

G AB dU A dU B = dF 2 + sin 2 F dVl 2 . (4) 

For the four dimensional spacetime M. we assume spherical symmetry and use the following ansatz for the metric 



g ab dx a dx b = e 2a ^ (-e 2 ^dt 2 + dr 2 ) + r 2 dw 2 . 



(5) 
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In addition we assume that the map U is spherically £-equi variant, that is 

F = F(t,r), fi = x«(w), (6) 

where Xii 1 - ) is a homogeneous harmonic polynomial of degree I. For this ansatz, the energy-momentum tensor @ 
does not depend on angles and thus can be consistently coupled to the spherically symmetric Einstein equations. We 
note in passing that a very similar idea of introducing the " angular momentum" into spherical collapse was put forward 

by Olabarrieta et al. Q in the context of critical phenomena. In terms of the mass function m(i, r) = 7j r (l — e~ 2a ) 

the Einstein equations take the following form (hereafter primes and dots denote partial derivatives with respect to r 
and t, respectively) 

m' = | r 2 e- 2a (V 2 + e^F 2 ) + J-^l sm 2 F , (7) 

m = Kr 2 e~ 2a FF', (8) 
2m o„ „,„ „, sin 2 _F 

2~ 



0> = ^e 2a -^+l)e 2a ^^, (9) 



where k = 4ttGX is a dimensionless parameter. The wave map equation ([3]) takes the form 

(e^F) ' - 1 (r 2 e^F')' + e^+ 2a €(£ + = . (10) 

For ^ = the above equations reduce to the Einstein-massless scalar field equations analyzed by us in [l[. For k = 
(no gravity) equations (7)- (9) are trivially solved by m — and (3 = 0, while equation (JTOJ) reduces to the flat space 
wave map equation. 



III. ITERATIVE SCHEME 



We assume that initial data are small, smooth, and compactly supported (the last assumption can be replaced by 
a suitable fall-off condition) 



F(0,r)=eg(r), F(0, r) = eh(r) . (11) 

We make the following perturbative expansion 

m(t,r) = m Q (t,r) + emi(t,r) + e 2 m 2 (t,r) + . . . , (12) 

P(t,r) = /3 (i,r)+£/3i(t,r)+e 2 /3 2 (i,r) + ..., (13) 

F(t,r) = F Q (t,r)+eF 1 (t,r)+e 2 F 2 (t,r)+e 3 F 3 (t,r) + .... (14) 



Substituting this expansion into the field equations and grouping terms with the same power of e we get the iterative 
scheme which can be solved recursively. 

We consider perturbations about Minkowski spacetime, so mo = /3o = Fq = 0. At the first order the metric 
functions m\ — f3\ — (this follows from regularity at r = 0), while F\ satisfies the flat space radial wave equation 
for the £-th spherical harmonic 

= , U (e) = d 2 - d 2 - 2 -d r + , (15) 

with initial data -Fi(0, r) = g{r), i*i(0, r) — h(r). The general everywhere regular solution of equation (|15p is given by 
a superposition of outgoing and ingoing waves 

F 1 (t,r)=F[ et (t,r)+Ff*(t,r), (16) 

where 

FV*(t r) = ly {2£ - k)l aik)[u) H*>(t r) = i V( i)^ ( 2£ - fc ) ! (i 7 ) 

1 l ' J r ^fc!(^-fc)!(«-u)^- fe ' 1 k ' J r ^ l j fc!(^ - fc)! (v - uf- k ' 1 ' 

and u — t — r 1 v = t + r are the retarded and advanced times, respectively (the superscript in round brackets denotes 
the fc-th derivative). Note that for compactly supported initial data the generating function a(x) can be chosen to 
have compact support as well (this condition determines a(x) uniquely). 
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At the second order Dm]^ = 0, hence F 2 = (because it has zero initial data), while the metric functions satisfy 
the following equations 



m> 2 = ^r 2 (F?+Ft 
m-2 = K r 2 F\F' x , 



£(£ + ! 



2 — F i I 



& = 



2m 2 £(1+1) 



F 2 



We temporarily postpone the discussion of this system and proceed now to the third order, where we have 

□ F — 1R F 4- R F A- R 1 F 1 2l(f+l)m 2 F 1 21(1+ l)Ff 

U (i)^3 - 2p 2 t x + P2-T1 + P 2 t x — 1 — . 



(18) 
(19) 
(20) 

(21) 



To solve this equation we use the Duhamel formula for the solution of the inhomogeneous wave equation O^F 
N(t, r) with zero initial data 



t t+r-T 



F(t,r) 



2r 



(It 



P Pi(fi)N(T,p)dp, 



(22) 



It — i — -r I 



where Pg(p) are Legendre polynomials of degree I and p = (r 2 + p 2 — (t — r) 2 )/2rp (note that — 1 < p < 1 within the 
integration range). Applying this formula to equation (|21l) , using null coordinates r\ — t — p and £ = r + p, and the 
abbreviation K(m, f3, F) = 2/3F + f3F + (3'F' - (2£(£ + l)/r 2 )(mF/r ~ F 3 /3), we obtain 



t+r t-r 



(23) 



where now p = (r 2 + (£ — t)(t — rf))/r(£ — rf). If the initial data (fTTj) vanish outside a ball of radius R, then for t > r + R 
we may drop the advanced part of F\ (t, r) and interchange the order of integration in (|23f to get 



00 t+r 



F 3 (t,r) 



8r 



dv / ^- V )PAp)K(m 2 (^ v ),f3 2 ^, v ) 1 Fr t (^v))^. 



(24) 



— 00 t— r 



In order to determine the late-time behavior of F%(t,r) we need to know the behavior of the source term K along 
the light cone for large values of r (the intersection of the integration range in (f24|) with the support of F[ et (t,r)). 
Having that, we shall expand the function K in ([24| in the inverse powers of p — (£ — rj)/2 and calculate the integrals 
using the following identity (see the appendix in 8} for the derivation) 



t+r 



; 2(n- 2 f- r e+1 (t-r]) n - e - 2 
~ (21+ 1)!! [(t - f]) 2 - r 2 ]™" 1 



F 



l+2-w l+3-rt 
2 ' 2 
£ + 3/2 



(25) 



where (n - 2)^ = (n - 2)(n - 3) • • • (n - I - 1). 

Now, we return to the analysis of the second-order equations (18)-(20). Substituting the outgoing solution (Tl"T|) into 
(TTSJ) and integrating, we get 



,, s t>R 



a^\t-p)) - £ E 



' + n)^(l + k - 1 - n) 2 ^- 1 -") (l 2 +£+(k- n)(n + 1)) 



9, 



l<fc<2£+2 0<n<fc-l 

aV+ 1+n - k \t - p) a^- n \t-p) 



k2 k (k-1- n)\n\ 



P 



dp, 



(26) 



where we used that rri2(t,r = 0) = 0, which follows from regularity of initial data at the origin and (|19[) . Here and 
in the following we use repeatedly the fact that a(x) = for \x\ > R, R being the radius of a ball on which the 
initial data (fTT|) are supported. To describe the behavior of ra 2 (£, r) along the lightcone it is convenient to use the 
null coordinate u — t — r instead of t, and rewrite (|26p as 



I \ r+u>R 

m2(u,r) = k 



2r 



aW( W )) : 



1(1 + \)(l 2 + I + 2) a^- 1 ) (u)a (£) («) 



O 



(27) 
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Next, using the gauge freedom to set @2{t,r — 0) = and integrating equation (|20|) . we get 

r oo 

ft(i,r) 2«y i J {a^ +1 \s)Y dsdp 



t-p 







(2g2 | 2g | , ^+l)(3£ 2 + 3^ + 2)^- 1 )(t-p) a W(t-p) | /l 

p 3 2 p 4 \p 5 



dp. (28) 



The first integral can be integrated by parts giving 



a I \ r+u>R „ 

P2(u,r) = 2k 



1 / («<«>«)= 



, (a^+ 1 )( S )) 2 

<& + / ? x <*« 

r — (s — it) 



(2£ 2 + 2^+1) 



(aW (a)) 2 + 1) (3^ 2 + 3£ + 2) a^ 1 ) («) W (a) 



(r — (s — u)) 3 



(r - (s - m)) 4 



+ -(29) 



IV. TAILS 



Now, we shall apply the method described above to compute the late-time asymptotics of solutions in the third-order 
approximation. Hereafter, it is convenient to define the following integrals (for non-negative integers m, n) 



(«)= !{s-uT (a {n) (s)j 



ds. 



(30) 



A. l = \ 



From we have 



r+u>R K 



r+ = >R -£ 



2ll{u)+ l -{2ll{u)-M»(u)) + o(\) 
2J 2 » + i(4Z2(u) - 5 (a») 2 ) + (1) 



2J 2 V)-:(a») 2 + Ol 3 



Substituting ([T7J) and (|3HI33|) into {2U we obtain 

hoc i+r 



F 3 (t,r) = y J dr, J d£ 



(Z-v) 2 



where 



tfxfo) - ^(^'(i?) + (2/ 2 2 (^) - 5/ 1 (r ? ))a / '(77) • 



(31) 
(32) 

(33) 



(34) 



(35) 



Performing the inner integral over £ in (|34[) with the help of the identity (|25p we get the asymptotic behavior which 
is valid for large retarded times u 



where 



(t 2 -r 2 ) 2 



+ 00 



kAi +0 



M = - I (a"(s)fa(s)ds. 



(36) 



(37) 



From p6p we obtain the late-time tails in both asymptotic regimes: F^(t,r) ~ nA%rt 4 at future timelike infinity 
(r = const, t — > oo) and (rF 3 )(v = oo, u) ~ nAi(2u)~ 2 at future null infinity (w = oo, u — > oo). 
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B. I > 2 



We give the detailed calculation only for £ — 2. In this case we have from (|2"9")1 



n , \ r+u>R H 

P2{u,r) = — 



2%(u) + -(2/ 2 (u) - 13J 2 °(u)) + -j(2J|(u) - 394(u) + 30 (a'(u)f) + 0[- 
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2i£(u) + -(47 3 1 (u) - 13 (a'») 2 ) + ^(6/ 3 2 (m) - 39/ 2 °(w) - 60a'(u)a'») + O 



2I°(u) - — (a"(u)f - %a'{u)a"{u) +o(\ 



(38) 
,(39) 
(40) 



Substituting (17]) and (|38H40p into we obtain 

+ 00 t-\-r 



F 3 (t,r) 



dr) / d£, 



-00 t— r 

4 



% ( J 3W a(3) ^ 



O 



where 



and 



(a"(rj)f + 2n [2 (a"^)) 3 - 3 (a^( V )) 2 a(rj) + + ^ l ( , _ . (/) , 

C/ 2 (r?) = BllWin) + (2/ 2 (, y ) - ll/ 2 °(ry))a( 3 )(r,) , 



(41) 
(42) 



Vafo) = -24J 3 (7 ? )a(^) + 36 J3 1 (77)^' (77) + (-117/°^) + 18/f (r?)) a"(r,) + (-781^) + 41* ( v ) + 60 (a'( V )f) a^( V ) 



(43) 

Performing the inner integral over £ in (|4ip with the help of the identity ([25]) we get the asymptotic behavior for large 
retarded times 



F 3 (t,r) 



{t 2 - r 2 f 



B 2 + 



where 



An 



128 
~15 



+ OO 



(a"(s)) 3 -3(a (3) (s)) 2 a(s) 



ds and Bo 



G4 
15 



{a"(s)fds. 



(44) 



(45) 



For the general I it is easy to see that the first nonzero contribution to the tail comes from the term with n — £ + 2 
in the identity (j2"B"|) which gives the following asymptotics 



F 3 (t,r) = 



(t 2 - r 2 )^ 1 



kAh + Bi 



(46) 



The formula (|46|) gives the first term in the asymptotic series approximation of the solution for late retarded times, 
that is for small e we have 



(t 2 



-\F(t,r)-e 3 F 3 (t )r )\=0(e 5 ). 



(47) 

r~ 

We have not attempted to derive a general formula for the coefficients Ai and Bg - the computation of these coefficients 
for each given £ is straightforward but as £ increases the algebra becomes tedious since it involves high-order expansions 
of the metric functions along the light cone. Anyway, it follows from (l46l) that the tail behaves as F 3 (t, r) ~ r £ £-( 2£ + 2 ) 
at future timelike infinity and as (rF 3 )(v = 00, it) ~ at future null infinity. 



Remark 1. For £ > 2 the tail (|46[) has two parts quantified by the coefficients nAg and Bg, respectively. The Ag-part 
comes from the gravitational self-interaction of the wave map and vanishes for k — 0. The /3^-part comes from the 
cubic nonlinearity of the wave map equation and is present without gravity as well. The case I = 1 is special in the 
sense that the -B^-part is absent in (f3"6")) since it is subdominant (decaying as t~ 5 ) with respect to the leading order 
term. 



Remark 2. It is instructive to compare the tail (|46|) with the tail for a test linear massless field propagating on a 
fixed stationary asymptotically flat background. According to the Price law 0, [13, [ll| the £-th multipole of this 
linear tail <j>(.(t,r) ~ r l+1 /(t 2 — r 2 ) i+2 for t — r — > 00. This decay is by one power faster than that in (|46|) . Of 
course, this difference is not very surprising as the tail studied here and Price's tail correspond to different physical 
situations, however we point it out as another example of the inapplicability of linearized theory in the study of 
radiative relaxation processes (see [13. [l3| for other examples). We shall discuss this issue in more detail elsewhere 
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V. NUMERICS 



In this section we compare the above analytic predictions with the results of numerical solutions of Einstein-wave 
map equations (|7] |10p for various initial data. The details of the numerical method were given in [l| for the case I = 0. 
The only difference for higher I is the boundary condition F{t, r) ~ r l for small r which guarantees regularity at the 
origin. The initial data were generated by the gaussian 



ea(x) = eexp (—x 2 ) 
for different values of e. For these initial data the formula (37) gives for I = 1 

G4 



Ax = —V3vr« 7.2769, 



and the formula (45) gives for £ = 2 



An 



10240 



3tt w -388.1, 



So 



2048 



3tt 



45.52. 



(48) 



(49) 



(50) 



81 ' 405 

In order to extract the parameters of the tails at timelike infinity we fit our numerical data with the formula 

F(t, r) = At 1 exp (B/t + C/t 2 ) . (51) 
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FIG. 1: Left panel: The log-log plot of F(t,r) for fixed r — 5. Fitting (15 1 we get power-law exponents 7 = —4.0196 {£ = 1), 
-6.0009 {£ = 2), -8.0049 {£ = 3), in agreement with the analytic prediction (46). Right panel: The log-log plot of F(t,r)/r e 
for fixed large advanced time v = t + r — 1200 as the function of retarded time u = t — r. The analogous fit to (51) yields the 
exponents -2.0036 {1 = 1), -3.0004 (I = 2), -4.0095 {£ = 3), in accordance with (46). In both panels k = 0.02 and e = 2.0 
{£ = !),£= 0.7 {£ = 2), e = 0.3 {£ = 3). 





A{£ = 


= 1) 




A(£ 


= 2) 


e 


theory 


numerics 


e 


theory 


numerics 


0.05 


9.096e-5 


9.051e-5 


0.05 


-0.07277 


-0.07274 


0.1 


7.277e-4 


7.289e-4 


0.1 


-0.58216 


-0.58476 


0.4 


0.04657 


0.04701 


0.2 


-4.65727 


-4.6841 


0.8 


0.37258 


0.37414 


0.4 


-37.2582 


-37.3778 


2.4 


10.0597 


10.0299 


0.65 


-159.875 


-160.441 


3.2 


23.8452 


16.0528 


0.7 


-199.681 


-189.377 


3.8 


39.9303 


19.6931 


0.75 


-245.598 


-189.792 



TABLE I: The comparison of analytic and numerical amplitudes of the tails at timelike infinity. Here k = 0.02 and r = 5. 
The third-order approximation reads A = e 3 nrAi for 1 = 1 and A = £ 3 r 2 {nA2 + B2) for £ = 2. 

The results and their confrontation with analytic predictions are summarized in Table 1 and Figures 1, 2, and 3. 
From this comparison we conclude that the third-order approximation is excellent for sufficiently small initial data. 
For large data approaching the black- hole threshold the third-order approximation breaks down - this is seen in Fig. 2 
as the deviation from the scaling A ~ e 3 and in Fig. 3 as the deviation from the linear dependence of A on n. 
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FIG. 2: The log-log plot of the amplitude of the tail at timelike infinity as a function of the amplitude of initial data (black 
dots) for fixed k = 0.02 and r — 5. The third-order approximation (dashed line) is excellent for small data, but it breaks down 
for large data lying near the threshold of black hole formation. 
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FIG. 3: The plot of the amplitude of the tail at timelike infinity as a function of the coupling constant k (black dots) for fixed 
e = 0.05 and r = 5. As c increases we leave the small-data regime and consequently the third-order approximation (dashed 
line) deteriorates. 



It should be emphasized that we get the same decay rates £~( 2£ + 2 ) (at timclikc infinity) and (at null infinity) 

for all subcritical evolutions, regardless of whether our third-order formula reproduces accurately the amplitude of 
the tail (for small data) or fails (for large data). 
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